Theoretical study of electronic Raman scattering of Borocarbide superconductors 
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The electronic Raman scattering of Borocarbide superconductors is studied based on the weak 
coupling theory with s + g-w&ve gap symmetry. The low energy behaviors and the relative peak 
positions can be naturally understood, while the explanation of the detailed shape of the B\ g peak 
seems to require a strong inelastic interaction not present in the weak coupling theory. 
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I. INTRODUCTION 



Electronic Raman scattering is a very useful probe 
in determining the symmetry of the order parameter 
of superconductors (SC).d They play especially impor- 
tant roles in the study of unconventional superconduc- 
tors owing to the strong dependence of Raman responses 
on the symmetry of superconducting order parameters. 
Along with the discovery of high temperature super- 
conductors, other unconventional superconductors have 
been actively investigated, for example, Borocarbide su- 
perconductors (BCSC)H (Y,Lu)Ni 2 B 2 C, and Ruthenate 
superconductorsfj Sr2Ru04. The symmetry of the order 
parameter of high temperature SC.and Ruthenate SC is 
believed to be d-waveo and p-waveO, respectively. How- 
ever, the symmetry of the order parameter of BCSC has 
not been determined unambiguously at present. Never- 
theless, there are growing experimental evidences for the 
strongly anisotropic nature of colder paramter of BCSCL 
The band structure calculationstm and T c versus 7 plotLl 
suggest strong electron-phonon interaction as a possible 
origin of BCSC. In general, the superconductivity me- 
diated by electron-phonon interaction has an isotropic 
s-wave superconducting gap (SG), and the nodes of gap 
on Fermi surface (FS) do not exist. 

On the other hand, the following experiments indi- 
cate the existence of the nodes of SG: (1) the power 
law .temperature dependence (T 3 ) of nuclear relaxation 
rateO 1 jT\ suggests the existence of line nodes of the su- 
perconducting gap on Fermi surface Jj Moreover, Hebel- 
Slichter peak is not observecB, which is also consistent 
with the existence of nodes. (2) The magnetic field (H) 
dependence of residual density of states in the vortex 
statcEjiid was observed to b&-f/H. According to the 
well-known result by VoloviktHl, the y/H indicates the 
existence of lines of gap nodes on FS. (3) V. Metlushko 
et al. reported an observation of the fourfold symmet- 
ric upper critical field.Ej Evidently, such a behavior is 
not compatible with the isotropic s-wave gap, and has 
been interpreted ia. terms of three dimensional d-wave 
superconductivity.!!-] (4) The direct measurement of SG 
by photoemission spectroscopy^ also demonstrates the 
existence of nodes of SG. (5) The recent experiment by E. 
Boaknin et al. on the thermal conductivity of LL1N12B2C 



shows that the magnetic field dependence of thermal con- 
ductivity at low temperature is very similar to that of. a 
heavy fermion SC with gap which has lines of nodesEH 

The above experimental facts (l)-(5) clearly corrob- 
orate a gapless superconductivity, and the simplest of 
which is d-wave superconductivity. However, the follow- 
ing experiments for the BCSC with Ni substituted by 
some impurity element as Pt pose serious challenges for 
the d-wave scenario: (6) The magnetic field H depep-, 
dence of residual density of states in the vortex statei9'til 
changes frona-j\/i? to H upon the impurity substitution. 
Caroli et alxH showed that the H dependence is linear 
for the isotropic s-wave SC. (7) The detailed plot of 
\n[C (T) / jT c ] versus T c /T also demonstrates the change 
of behavior from T 3 to the activated one£3 (8) The pho- 
toemission experiment shows the opening of gap over the 
whole FS for impurity substituted compound by compar- 
ing the shapes of the spectra of pure and dirty BCSC. 

The experiments for the dirty BCSC (6)- (8) imply 
the gap symmetry cannot be a simple d-wave because 
the introduction of impurities does not change the gap 
anisotropy owing to the d-wave symmetrytS, and only 
overall amplitudes decrease with the impurity concentra- 
tions. In connection with the experiments for the diry 
BCSC we note that the superconducting temperature of 
conventional s-wave SC is insensitive to the introduction 
of non-magnetic impurities, which is the fact known as 
Anderson theorernEEl, while, the unconventional SC de- 
pend on them sensitively.EElEJ 

Finally there is no experimental evidence of broken 
tetragonal symmetry, thus we are forced to rule out pair- 
ing interactions which violate tetragonal symmetry. 

Synthesizing all of the above observations, (l)-(5), (6)- 
(8), and the unbrokent tetragonal symmetry, we are led 
to a proposal that the gap symmetry of BCSC should 
be a mixed one which respects the tetragonal symme- 
try. Obviously, the simplest possible candidate is the 
s + g-w&ve gap symmetry. The gap function of the form 
A(^) = |Ao cos(20)| is also a viable choice, but the nodes 
of the gap of this kind are accidental and are not enforced 
by symmetry. From the symmetry viewpoint, s + g-wave 
gap is more natural. .— . 

Yang et al. studiedB BCSC by means of electronic 
Raman scattering. They observed 2A-like peaks in 
Aig,Bi g , and Bi g geometry. The peak position of B\ g 
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geometry is found to be larger than the others. Also, 
the peaks of the Ai g ,B 2g symmetry are much stronger 
than the peak of B\ g symmetry. There exists scattering 
strength below 2A gap, which increases with frequency 
linearly in B\ g ,Bi g geometry. [See Fig. 5] 

We intend to understand the qualitative features of the 
above Raman experiments based on the simplest model 
with s + g- wave gap symmetry. We adopt the weak cou- 
pling theory with separable pairing interaction. Then, 
even if the physical origin of pairing interaction is un- 
known, we can deduce many physical consequences which 
mainly depend on the pairing symmetry not on the de- 
tails of interactions. Strictly speaking, the weak cou- 
pling theory is not expected to valid on a quantitative 
level, and indeed the band theory calculation predicts the 
electron-phonon interaction to lie in the strong coupling 
regime. Nevertheless, we believe that the qualititave fea- 
tures might be understood with the simple weak coupling 
model. 

We find that the relative peak positions and the nature 
of scattering strengh below 2A gap can be understood 
naturally, while the weak and broad B\ g peak features 
cannot be explained in a simple weak coupling theory 
with s + g-w&ve gap symmetry. Presumably, in B\ g ge- 
ometry, other inelastic scattering mechanisms in particle- 
hole channel become important. The B\ g phonon may 
provide one of such scattering mechanisms. 

This paper is organized as follows: In Sec. II, we briefly 
summarize basic formalisms of electronic Raman scatter- 
ing. In Sec. Ill, the explict expressions of Raman inten- 
sity of s+g-w&ve SC are derived. The physical results are 
presented in Sec. IV, and we conclude this paper with 
summary and discussion in Sec. V. 



II. FORMALISM 

The general formalisms for the electronic Raman scat- 
tering of superconductorSj-are well expounded in the pa- 
pers by Klein and Diekeifil, Monien and ZawadowskicJ, 
and Devereaux and EinzelcJ. Here we shall briefly sum- 
marize the main results of the above papers necessary for 
our discussions. 

The electronic Raman cross section is proportional to 
the dynamical structure factor S(LO,q): 



S(w,q) = [l+n fl (w)] 



- Im xfp(w, q) 



(1) 



where tib(ui) is the Bose distribution function and the su- 
perscript R of xfp denotes the retarded correlation func- 
tion. The effective density-density correlation function 
Xpp (in imaginary time) is defined by 

Xpp(r-r', q ) = -(T T p q (r)p_ q (r')). (2) 
pq is the effective density operator 



7 P c P+q/2^ c p-qA<r- 



(3) 



For non-resonant electronic Raman scatterings, the coef- 
ficient 7p is given by 



7p 



(4) 



where e ,e F is the polarization vector of the incoming 
and outgoing photons, respectively, ep is the energy dis- 
persion of material, a, (3 label the coordinates perpen- 
dicular to the photon momentum. 

Light couples with the charge fluctuations in the ma- 
terial, and in this case the long-range Coulomb inter- 
action should be summed in random phase approxima- 
tion (RPA) . The RPA summation implements the screen- 
ing effect. Let us define irreducible correlation functions 
which do not include the RPA type diagrams as subdia- 
grams. 



7T 77 (q,T-T') 
K0~/(<l,T- t') 

7r 00 (q, t — t') 



where 



Pel 



= -( T rPq(T)p-q(r'))i rr , (5) 

= -(7Vpq(r)p-q(r'))i rr , (6) 

= -(T r/9 q(T)p_q(r')>i rr , (7) 

= -(T T/9 q(r)p_q(r')) irr! (8) 

t 



E p .<t 1 c p+q/2 l( T c p-qA<r 



is an ordinary 



density operator. Then the correlation function \pp can 
be written as 



Xpp 



7r 70 7To 7 



Too 



7r 7 o7To 7 



Xpp j 



(9) 



[ oo 



where Xpp IS the density-density correlation function 
which is negligible in the low q limit. In the low q limit, 
only the first term in the bracket of Eq. (^) needs to be 
considered. 



qi™ Xpp 0'w,q) 



'77 



(iu),q) 



7r 7 o(iw,q)7To 7 (zw,q) 



7r 00 (iw,q) 

(10) 

The second term of the right hand side of Eq. (110) is 
usually referred to as the screening correction. 

The vertex factor 7p for the lattice with tetragonal 
symmetry in various scattering geometries are given by 



lA lg {4>) 
lB lg {4>) 
7B 2g (0) 



1 + 7a cos 40, 
7Bj cos 24>, 
7s 2 sin 2(f), 



(11) 



There can be terms with higher harmonics such as 
cos(4A0) 1 cos((4A - 2)0),sin((4A - 2)0) with N > 2 
m Eq. (P). We will ignore them, but in the detailed 
fitting of experimental data, the higher harmonics terms 



are necesarycij The factor of 1 for the A\ g geometry case 
of Eq. ( p"l| ) is cancelled by the screening correction, there- 
fore, effectively 7Ai s (0) - * J A cos 40. 

The coefficients 7_a, 7Bi , 7b 2 are determined by the 
shape of energy band. If we take a simple band structure 
with tetragonal symmetry [ a is lattice constant] 

ep = — 2t[cos(p x a) + cos(p y a)] + At' cos(p x a) cos(p y a) 
- 2t z cos( Pz a). (12) 
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then 7s ls oc t , jB 2g x t' an d 1a oc const. t + const. t'. In 
our work, we will simply put 1a — 1Bi = 1b 2 = !• Thus, 
there are some reservations in comparing the relative in- 
tensities of theoretical results with those of experimental 
ones. 



III. THE s + 9 - WAVE SUPERCONDUCTOR 



The gap function of (s + g)-wave superconductor is 



A(<f>) = A s + A fl cos(4<£), 0<A s <A ff , 



(13) 



where <f> is the azimuthal angle in spherical polar coordi- 
nate. The weak dependence on the radial momentum is 
neglected (k = kp). Since the amplitude of <?-wave part 
is larger than s-wave part the lines of nodes exist on the 
FS. In the opposite case A s > A g the gap function does 
not have nodes, and it behaves like gapful anisotropic 
s-wave SC. The case < A s < A g and < A g < A s 
is physically relevant for pure BCSC and dirty BCSC, 
respectively. We will use the Nambu formalism for the 
computation of correlation functions in superconducting 
state. The one-particle Green function in weak coupling 
approximation is given by 



G(ie,p) = 



i£T + £pT 3 + ApT"! 



A 2 



(14) 



where Tj's are the Pauli matrices in Nambu space. The 
wave function renormalization Z(ie) is neglected, and the 
frequency dependence of interaction is encoded in the 
cutoff w c [SeeEq. ©] of pairing potential in momentum 
space. In Nambu notation, the correlation functions 7r 77 ' 
can be expressed as [ p — [it, p), q = (iuj, q)] 



7rr r 



(«) = Jf E Tr [^p+Mp) i'p+qGip + q)] , (15) 



where N is the number of lattice sites. r PjP + q is the 
renormalized vertex function. We will evaluate the vertex 
function in ladder approximation. The unrenormalized 
(bare) vertex function is 

fp, P + 9 = 7p = T37p- (16) 
The gap equation of (spin-singlet) s + g-wave SC reads 

AM 



N 



P' 



2Er 



where £*pp' is the pairing potential in momentum space, 

and E p , = v /ep7TA p 7. 

The gap equation Eq. ([l7]) can be solved in a closed 
form for separable pairing interactions. Given the specific 
form of gap Eq. (E3T), the separable pairing interaction 



which respects the tetragonal symmetry should be chosen 
as 

£>pp' = [D s + D g cos A(j> cos 4</>' 
+ £> S9 (cos40 + cos40')] 
x e(w c -|ep|)e(w c -|$p/|), (18) 

where Q(x) is the step function. The numerical values 
of the coupling constants D s , D sg , D g are determined to 
be consistent with the values of A s , A g , and the super- 
conducting transition temperature T c . For details, see 
Appendix C. Since there is no dependence on the polar 
angle 9 in Eq. (|l5|), the momentum integration can be 
written as 



1 f 1 f 2n 



(19) 



where Np is the density of states per spin. Note that the 
spin sum is implicit in the trace over Nambu space. For 
later references, let us define the dimensionless coupling 
constants. 



A s = -N F D S , X g = —NpD„, X sg = —N F D 



sg ■ 



(20) 

The technical analysis of Eq. (|l^) crucially depends on 
whether the bare or the renormalized vertex is used. 



A. Analysis for the bare vertex 

Henceforth, we will consider only zero temperature 
case for simplicity. For the bare vertex, the integrals 
over frequency and |p| of Eq. (na) can be done explicitly. 
Then the imaginary parts of 7r?E/ can be expressed as a 



single integral over azimuthal angle 



Inmi 



-Re 



#7(</>) 7 '(0)A 2 (<£) 



o 2tt ,/uj 2 - 4A 2 (0) 



(21) 



At finite temperature the factor tanh(w/4T) should be 
multiplied to the above result. 

At this point we note that the screening correction is 
absent for B\ g and £?2 ff symmetry. The screening correc- 
tion can re-expressed as 







2k 



2tt 



[7W>)-l]A 2 (^)/(cos4^, 



(22) 



where / is a certain complex function which can be 
expanded in a power series with respect to argument. 
Clearly the integrals of Eq. (^2|) vanish for the vertex 
j(4>) = cos2</> or 7(0) = sin 20, since they are orthogonal 
to cos AN (f) with N > 1. Thus, the screening correc- 
tions for Bi g and B2 g symmetry vanish. The A± g vertex 
cos(40) has a finite overlap with the rest of integrand, 
and there is non-zero screening correction for A\ g sym- 
metry. The integrals Eq. (^lj) cannot be done in closed 
form, and the numerical integrations are required. The 
results are presented in Sec. IV A. 
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B. Analysis for the renormalized vertex 

The vertex correction will be evaluated within ladder 
approximation. Wc will include the corrections coming 
from the pairing interaction only. In principle, other in- 
teractions may well contribute to the vertex correction. 
However, we do not have enough experimental informa- 
tions to identify the specific interactions responsible for 
scatterings in particle-hole channel. Evidently this ap- 
proximation should be improved, and in fact, from the 
comparison with the Raman experiment we expect the 
existence of strong inelastic scatterings with B\ g sym- 
metry in particle-hole channel. Here we choose to work 
in the simplest approximation of including the pairing 
interaction only. 

The summation of ladder diagrams is equivalent to 
solving an integral equation of the following type: 

f 7 (p + q/2, p - q/2, iu) = 7 (p + q/2, p - q/2) 

E [(-l)^p P ']T3G(p'- q /2,ie') 



ie',p' 



xf 7 (p' + |, p' - G(p' + ^,ie' + ioj) r 3 .(23) 

The retardation effect in weak coupling approximation 
is reflected in the cutoff of pairing interaction in mo- 
mentum space, and then the frequency channel can be 
treated as instantaneous. For such interaction, the ver- 
tex correction T becomes independent of incoming fre- 
quency ie, which is an extremely crucial simplification. 
In fact, for a fully self-consistent treatment, the renor- 
malized Green function G with (Z(ie) ^ 1) should be 
used, which incoporate the self-energy correction due to 
the the scatterings in particle-hole channel. However, if 
the underlying system is Fermi liquid, Z(ie) correction is 
not expected to introduce qualitative changes. 

The momentum transfer f±, jeai3-.be neglected for SC 
with large penetraton depth.HE-IEa In the limit q = 0, it 
can be shown tha1n3 

f 7 (p,iw) = r 27 (p,icj)r 2 + r 37 (p,iw)T 3 . (24) 

Then Eq. simplifies to \p' = {ie',p')] 

r 27 (p,Kj) = — 2J [-T2-,(p',iu)A(p',iuj) 

P',W 

+ r^(p',iw)B(j/,iu)] x {-D pp ,} 

rp p 

r 37 (p,wj) = 7p + ^ 2J Y-^{v',iu)C(p',iu) 
P',«' 



B(p',iL0) 



ie'(ie' + iio) - A^, + 
iAyy iio 

Wf^W^+W^)' (26) 



The symmetry structure of pairing potential Eq. (|1S| ) 
and the separable kernel form of Eq. ( p5[ ) imply that the 
vertices for Bi g and B 2g symmetry are not renormalized, 
namely: 

r 27 (p,iw) = 0, r 37 (p,iw) = 7 p , for B lg ,B 2g . (27) 

If the additional scatterings were present in particle-hole 
channel with _B lff and/or B 2g symmetry the vertices for 
Bi g and B 2g symmetry would be renormalized. 

The non-trivial correction influences only the A\ g ver- 
tex without the additional interactions mentioned above. 
The vertex correction for the A\ g can be calculated ex- 
plicitly. The details of the calculations can be found in 
Appendix B. 



IV. RESULTS 
A. Results with the bare vertex 

Even though the analytic solutions in closed forms are 
not available, the asymptotic behavior at low frequency 
and the peak positions of correlation functions can be 
understood without numerical integrations. [For a pure 
d-wave SC, the integral can be done exactly in terms 
of elliptic integrals.E3] The detailed analytic calculations 
of low energy behaviors of Raman susceptibilities 7r 77 
can be found in Appendix A. Numerical integrations are 
required for the solutions valid over the entire frequency 
range. The numerical results are presented in Fig. [j] and 
Fig. [|. The result shown in Fig. [j] is for the case A ff is 
slightly larger than A s , while Fig. g shows the result for 
the case A g is much larger than A s . 

The Bi g spectrum is characterized by a sharp peak at 
lu = 2(A S + A g ), and a linear frequency dependence at 
low frequenc y. T he low frequency behavior is described 



by [See Eq. ( J A 1 3| ) 1 (1 - A s /A ff )w + const.w 3 . Note that 



the slope of linear term is larger for smaller A s / A g . 

The A\ g and B 2g spectra are characterized by a peak 
at lu = 2(A g — A s ), and a linear frequency dependence 
at low frequency. Without the screening correction the 
Ai g spectrum would have peaks at both to — 2(A S + 
A g ) and oj = 2(A g — A s ). [See discussions below.] The 
peak at oj = 2(A S + A g ) of A\ g spectrum is cancelled 
by screening correction, however, a small hump feature 
+ T 27 (p ,iio)B(p ,iuj) x [— -Dpp'], (25) remains at the frequency. B 2g spectrum does not show 



where 

A(p',iuj) = 



((z e ') 2 -^)((^ + H 2 -^)' 



noticeable structutre at to = 2(A S + A g ). 

The peak positions of the above spectra are determined 
by the following three conditions. (A) The angular in- 
tegrals Eq. ( El] , |22] ) can be re-expressed in terms of 
z = cos4</>, then the Jacobian of the change of variable 
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— z 2 appears in the integrand. Then, large con- 
tributions to the integrals stem from the angular region 
where z = cos 40 ~ ±1. This feature is very similar to 
that of a van Hove singularity. (B) Large contributions 
to the integrals also come from the angular region where 
lu 2 - 4A 2 (0) - 0. [ See Eq. © ] The frequencies which 
satisfy the conditions (A) and (B) are lu = 2(A g ± A s ), 
and the peaks are expected at those frequencies. (C) The 
conditions (A) and (B) predict peaks at both frequencies 
lu = 2(A g + A s ),2(A s — A s ) irrespective of the scattering 
symmetry. The scatttering symmetry is encoded in the 
vertex factor 7(0). 

For B lg symmetry, 7 2 (0) = cos 2 (20) = 1+c ° s40 . The 
vertex ( squared ) vanishes at cos 40 = — 1 , and this factor 
suppresses a peak that would appear at lu = 2(A g — 
A s ). [ Recall A(0) = A s + A g cos40 ] Therefore, for 
Big symmetry, only a single peak at lu = 2(A g + A s ) is 
expected, which agrees with the numerical result. 

The vertex factor of Bi g is 7 2 (0) = sin 2 (20) = 
i-cos 40 ^ Following t ne same arguments as the Bi g case, 
we may expect a single peak at lu — 2(A g — A s ), which 
also agrees with the numerical result. The screening cor- 
rections do not alter the above results because the cor- 
rections simply vanish for B\ g ,B2 g spectra as discussed 
in Sec. Ill A. ' 

For the A lg case, we have to take the screening correc- 
tion into account. Now notice that in the angular region 
cos 40 ~ 1, [ SeeEq. (|ll|) ] 

K~tA 1 lA 1 ~ 7>7<4i0 ~ "00 ■ (28) 

Then, the contributions to the Raman susceptibility, 

Xia 1 ia 1 = k 1Ai ~ 1Ai - (TT JAi o) 2 /n n, (29) 

coming from the angular region cos 40 ~ 1 becomes small 
due to the cancellation between the first and the second 
term. The above cancellations also occur for the con- 
tributions from the angular region cos 40 ~ —1. Then, 
subleading contributions determine the peaks. Numerical 
results show that the A lg peak appear at lu — 2(A 9 — A s ). 



B. Results with the renormalized vertex 

As discussed in Sec. Ill B, the vertex correction affects 
only Ai g spectrum. The vertex correction depends on 
the relative magnitude of A s and A g . The results are 
presented in Fig. |^ and Fig. ||. It turns out that the 
vertex correction is more important for larger A g /A s . In 
case of Fig. || [A g /A s — 1.28], the vertex correction is 
almost negligible over the entire frequency range. On the 
other hand, in case of Fig. [| [A g /A s = 5.8], the vertex 
correction suppresses the peak at lu — 2(A g — A s ) which 
was present in the spectrum with the bare vertex. At low 
frequency, the vertex correction is negligible as in the case 
of Fig. ||, while at frequency higher than 2(A 9 + A s ), it 
slightly enhances the spectrum. More discussions will 
follow in the next section. 



Intensity 
Arb. Units. 




1 2 3 4 5 6 7 



Frequency(mev) 

FIG. 1: Raman Intensities with the bare vertex. A s = 
1.72mev, A s = 1.34mev, and T c = 15. 9K. 



Intensity 
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1 2 3 4 5 6 7 

Frequency(mev) 



FIG. 2: Raman Intensities with the bare vertex. A g = 
2.5mev, A s = 0.43mev, and T c = 14.2K. 



Intensity 



Arb. Units. Alg Spectra 




1 2 3 4 5 6 7 

Frequency(mev) 



FIG. 3: Raman intensity with/ without vertex correction 
for Ai g spectra. A s = 1.34mev, A 9 = 1.72mev, and T c = 
15. 9K. X s = 0.13, A 9 = 0.22, X sg = 0.15. The solid line is the 
spectra with the vertex corretion, and the dotted line is the 
bare vertex spectra. 
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V. DISCUSSIONS AND SUMMARY 

Comparing the theoretical results, especially Fig. || 
, with the experimental data Fig. [5] (b), we find that 
the relative order of peak positions and the low fre- 
quency behaviors coincide with those of experimental 
data. In our model, the peak of B lg spectrum is located 
at uj = 2(A S + A g ), while that of B 2g , Ai g is located at 



2(A fi 



A„ 



Then A g and A s 



can be extracted 



from the experimental data, and we obtain A g ~ 2.5mev 
and A s ~ 0.3mev. Thus, it turns out that the s-wave 
component of gap is much smaller than the g-wave com- 
ponent. In fact, Fig. || [A g /A s = 5.8] compares more 
favorably with the experimental data Fig. |E| than Fig. [j] 
[A s /A s = 1.28] . 

We have to note that the linear frequency dependence 
Bi g spectrum at low frequency observed in the Raman 
experiment rules out the pure d-wave scenario, which 
predicts lo 3 dependence. [See Appendix A] The low fre- 
quency behaviors are robust againt the effect of vertex 
correction and provide an unambiguous evidence for the 
s + g-wave scenario. 

But the theoretical results pose some drawbacks. Most 
significantly, the shape of B\ g peak looks very different 
from the experimental one. Theoretical results predict 
a strong sharp peak at lo = 2(A g + A s ). Experimen- 
tally, the Bi g peak is very broad and weak compared to 
Ai g ,B 2g peaks. 

However, there exist some caveats in comparing the 
theoretical results directly with the experimental data. 
Recall that we have included only the pairing interaction 
in particle-hole channel in computing the vertex correc- 
tion. Owing to the specific symmetry of pairing interac- 
tion, the Bi g and B 2g channels were immune from the 
vertex correction. In general, we may well include the 
interactions in Bi g and B 2g symmetry in particle-hole 
channel which are compatible with tetragonal symmetry 
such as V cos 2<f) cos 2cf>', W s'm2(f>sm2(j)' . If such interac- 
tions were present the Bi g and /or B 2g channels surely 
would be influenced by vertex corrections. 

Regarding the discrepancies on the shape of Bi g peak 
structure between the theory and the experiment, we pre- 
sume that there exists a very strong inelastic scattering 
acting in Bi g channel which suppresses and broadens the 
Big peak. The introduction of inelastic scattering acting 
in Big channel would not modify the low frequency be- 
haviors since they are mostly dictated by the symmetry 
of superconducting order parameter. 

Borrowing the arguments in Sec. IV A, we may ex- 
pect that the effect of inelastic scattering in Bi g channel 
would be strongest near lo — 2(A g + A s ). Then, the weak 
peak structure at lo — 2(A 9 + A s ) of renormalized Ai g 
spectra [See Fig. 4] would be suppressed, and a peak 
would show up again at lo = 2(A 9 — A s ) in agreement 
with the experiment. 

The physical origin of the hypothetical inelastic scat- 
tering acting in Bi g channel is not clear, but the Raman 
experiment [ Fig. 5 (a) ] strongly suggests that the scat- 



tering between Bi g phonons and electrons might be a 
possible candidate. 

We also note that the inclusion of the detailed mo- 
mentum anisotropy of SG compatible with the tetragonal 
symmetry would give the results which will agree better 
with the experiments. 

In summary, we have studied the electronic Raman 
scattering of Borocarbide superconductors based on the 
weak coupling theory with s +g- wave gap symmetry. The 
low frequency behaviors and the relative peak positions 
can be naturally understood within a simple one-loop cal- 
culation. However, there are discrepancies in the detailed 
shape of peak structures between theoretical predictions 
and experimental results. The inclusion of the strong in- 
elastic scattering between Bi g phonon and electron may 
resolve the discrepancies. 
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APPENDIX A: LOW FREQUENCY BEHAVIORS 
OF CORRELATION FUNCTIONS 

At low frequency, the angular integral of Eq. ( |2l|) is 
dominated by the contributions from the angular region 
where A(<fr) ~ 0. In the region, the gap function A((f>) 
can be linearized with respect to the angle c/)q defined by 
A(0 O ) = 0. 

Let us first consider a pure d-wave case A(<f>) = 
Ad cos 2(f>. If we take into account the contributions 
from the first quadrant only, the gap function vanishes 
at (j) = 7r/4 . Writing <fi = tt/A + x, \x\ <C 1, the integral 
of Eq. (21) can be reduced to 



rim * f** d X l \cp=l+x)[AA^] (M) 



'•' 1 2ir ^2 ~ 16A 
The vertex factors become [see Eq. (|ll 

-1, 

-2x, 



2„2 



7B 2s (^ +x) ~ 1. 



Using scaling x = j^-y, we obtain for lo <g; 2A^ 



(A2) 
(A3) 
(A4) 
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-Dig 



-Ini7Tr 
4m7r s,„ ( w ) 



a; 



(A5) 



where the superscript (D) indicates the symmetry of su- 
perconducting order parameter. The results Eq. QAa) . 
coincide with those obtained by Devereaux and Einzcl.t£l 
Next let us consider the s + <?-wave case A(0) = A s + 
Upon linerization, one can write 



A g cos 



A(<f> — (f>o + x) ~ —ax, a = 4A ff sin4<^o, 
where the angle 0o is defined by the relation 

A s + A s cos 40n = 0, A s 
Then, the angular integral can be expressed as: 
N F dx 7 2 (0 = 0o + x)[a 2 x 2 ] 



,A g >0. 



Im7r 



R 



CO 



2tt 



Aa 2 x 2 



The vertex factors become 



7a 1s 
7iL 



6 + x) 



t>0 + x) 



cos" 



Ao = (A s /A s ) 2 



7b 2o (0o + x) 



1 - A s /A 9 
2 

[4 cos40 o ]x 2 
1 + A s /A g 



[2 sin40 o ]x 



(A6) 
(A7) 

(A8) 

(A9) 

(A10) 
(All) 



In Eq. ( A10 ) the second term is odd in x, thus it gives the 
vanishing resu lt up on integration. The first and the third 
term of Eq. ( |A1C| ) dominate the low energy behavior. 
Carrying out the integrals by rescaling as in d-wave case 
we get 



Im7r 



Im7r 



R(SG) 
Ai g 

R{SG) 



is 



j 2 Ai (A s /A g ) 2 co, 



1 - A s /A 9 



lm7T B { 2 s g G] 'M 



const. o; , 

2 1 + A a /A fl 1 
7b, [ s ] w ' 



(A12) 

(A13) 
(A14) 



where the superscript (5G) indicates the symmetry of 
order parameter. The coefficients of vertex factors 
7Aij 7bi> 7b 2 are reinstated for clarity. The above re- 
sults differ from the pure d-wave case by the w-linear 
component of B\ g spectra. We also note the dependence 
on the gap ratio < A s /A g < 1. For the small gap ra- 
tio, the w-linear component of A\ g spectra is suppressed, 
while that of B\ g spectra is enhanced. This agrees well 
with the numerical result. [See Fig. [l] and Fig. ^ 



APPENDIX B: CALCULATION OF VERTEX 
CORRECTION 



First, write 

r2,3(iw,p) = r? 3 (io;) + Tf 3 (zo;)cos' 



(Bl) 



Define two-component column vectors, 



T 2 (iLu) 



T s 2 (iuj) 

rf(w) 



> r 3 (zw) 



T s 3 (iuj) 

rg(i W ) 



(B2) 



Then, the integral equation for the vertex correction Eq. 
(Pa) can be recast in a matrix form. 



(B3) 



where h is the 2x2 identity matrix, and A, B, C are 2x2 
matrices to be defined below. 7 is a two-component col- 
umn vector defined by 



h + A -B 




"r 2 " 




' " 


-B I 2 -C 




. r 3_ 




1 



7 
7 



, for 7 = 1, 

, for 7 = cos 40. 



(B4) 



The matrix equation Eq. (B3) can be solved, yielding 
the desired solution: 



27 



3 7 



{h-G)B- 1 {h+A)-B\ 7, 
{h-C)-B{h+A)- 1 B\ 7. (B5) 
The 2x2 matrices A,B,C are defined as follows: 

A= \ ai H, B= H 1 I 2 }, C= \ Cl C2 1. (B6) 
a 3 a 4 J [ 63 64 J ' [ c 3 c 4 

The matrix elements are given by \p — (ie, p)] 

((j>)A(p,iu)), 



a 2 (iuj 
as(iuj 
04(10; 
bi(iu> 
b 2 {iuj 
6 3 (io; 
04(20; 
C\(iu) 



ie,p 



ie,p 



T 

N 

T 

N 

T 
N 

T 
N 



ie,p 
ie,p 
ie,p 
ie,p 
ie,p 

]vE 7 9 

ie,p 



0)(cos 40)A(p, io;), 

0)(cos40)A(p, ioj), 

4>) B (p, iuj ), 
0)(cos 4<f))B(p, ioj), 
(j))B(p,iuj), 
0)(cos 4(p)B(p, iu), 
4>)C{p, iu), 
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ie,p 

ca(iw) = — 23 I g (<j))C(p,iuj), 

d(ico) = I B (4>)(ooa44>)C(p,iu), (B7) 



A s (A g /i + A S9 / ) + A g (A g / 2 + Xsgh), 



dtp 



— (cosx) — ^ . 



(CI) 



Solving Eq. (CI) at T = T c , we obtain a relation between 
T c and coupling constants. 



where 



I s (4>) = X s + A sg cos 40, I g (<f>) = X g cos 40 + A. sg , 
A s = -iV>I>„ A ff = -N F D gi X sg = -N F D sg . (B8) 



The vertex correction Eq. (B5) together with the Eq. 
(Bl) should be substituted into the integral Eq. (15). 
Note that the solutions Eq. (B5) depend only on the ex- 
ternal frequency iu, so that just an analytic continuation 
iuj — ► uj + iS without additional integraions is sufficient 
to obtain the final answer. 



APPENDIX C: CHOICE OF COUPLING 
CONSTANTS 



lll[L1 35^] = <i±«-[^- A -> 



2 ' ^% 



1/2 



A s A g - X 2 sg 

(C2) 

For a pure g-wetve case, the well-known result T c — 
1.135 x Lu c e~ 2 / Xg is reproduced. Fixing the cutoff uj c — 
lOOmev, we have chosen the follo win g two sets of the 
coupling constants by solving Eq. (CI, C2). 



(I)Aa 



0.13, A 



0.22, A S3 = 0.15, 



1.34mcv, Ag = 1.72mev, T c = 15.9K.(C3) 



The gap equation Eq. ( |l7|) can be written as a set of 
coupled equations. Notations are defined in Eq. (B£). 



(II)A S = 0.05, Ag = 0.42, A sg 



0.06, 



A s = 0.43mev, A g = 2.5mev, T c = 14.2K. (C4) 



A s (A s / + Xsgh) + Ag(A s + A S g/ 2 ), 
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FIG. 4: Raman Intensity with vertex correction for A\ g 
spectra. A s = 0.43mev, A 9 = 2.5mev, and T c — 14. 2K. \ s — 
0.05, A 9 = 0.42, X sg = 0.06. The solid line is the spectra with 
the vertex corretion, and the dotted line is the bare vertex 
spectra. 




FIG. 5: Raman spectra-jf YNi 2 B 2 C with T c = 15. 3K from 
Yang et a/.'s experiments. The strong peak of (a) near 200 
cm -1 is due to B lg phonon. 



